Pathway interference in a loop array of three coupled microresonators 
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A system of three coupled toroidal microresonators arranged in a loop configuration is studied. 
This setup allows light entering the resonator setup from a tapered fiber to evolve along a variety 
of different pathways before leaving again through the fiber. In particular, the loop configuration of 
the resonators allows for an evolution which we term roundtrip process, in which the light evolves 
from one resonator sequentially through all others back to the initial one. This process renders 
the optical properties of the system sensitive to the phases of all coupling and scattering constants 
in the system. We analyze the transmission and reflection spectra, and interpret them in terms 
of interference between the various possible evolution pathways through the resonator system. In 
particular, we focus on the phase dependence of the optical properties. Finally, we discuss possible 
applications for this phase sensitivity induced by the roundtrip process, such as the measurement of 
the position of a nanoparticle close to one of the resonators, and the measurement of changes in the 
refractive index between two resonators. Our analytical results for the applications are supported by 
proof-of- principle calculations based on finite-difference-time-domain solution of Maxwell's equations 
in two dimensions on a grid. 

PACS numbers: 42.60.Da,42.82.Et,42.25.Hz,42.50.Ct 
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I. INTRODUCTION 



In recent years, optical microresonators have received 
considerable attention, since they offer a wide range of 
applications such as strong-coupling cavity quantum elec- 
trodynamics, the modification of spontaneous emission, 

optical communication, or as sources of light [Eg. a 

particular promising example combining several of these 
ideas is the goal of establishing quantum networks [!, 0] ■ 
By now, a large variety of implementations has been 
achieved 0, [Bl-lll| . Naturally, also the extension to more 
than one cavity has been suggested, for example, as 
chains of coupled ring or disc resonators [IH, [l3j] , of de- 
fects in photonic crystal hosts [1J, [l5[, of coupled res- 
onator spheres (lq - [l8j . or of coupled square resonators 
on a grid [l9j . 

A particular variant of coupled cavities involves two- 
dimensional arrays of microcavities, which can be used, 
e.g., to form photonic molecules [20. l2lj. or optical fil- 
ters [13, l22l - |24j . Light entering such an array can take 
a number of different pathways inside the cavities before 
leaving the coupled system, and the interference between 
these different pathways determines the optical proper- 
ties of the resonator system. This in a certain sense can 
be seen in analogy to an atom with multiple energy lev- 
els connected by several driving laser fields. Also the 
atom can evolve via different pathways. But a particu- 
larly interesting case arises if the laser fields are applied 
to the atom in a so-called closed-loop configuration [25r - 
Hlj]. This means that the laser fields are applied such 
that the atom can evolve in a non-trivial loop pathway 
from one initial state through the level scheme back to 
the initial state, e.g., |1) — > |2) — > |3) — > |1) with atomic 
states \i) (i £ {1,2,3}). The loop structure induces rich 
possibilities for interference between the different path- 
ways, and at the same time renders the optical properties 



of the atoms sensitive to the phase of the applied driv- 
ing fields. This prompts the question whether similar 
interference effects and phase-sensitivity could also arise 
in arrays of microresonators. The couplings between the 
resonators are mediated via their evanescent fields, and 
the corresponding coupling constants are in general com- 
plex. Similarly, the scattering inside a given cavity is 
characterized by a complex scattering constant. Thus it 
is not surprising that even for simple systems such as two 
coupled cavities the phase of the couplings can strongly 
influence the output fluxes. However, there are also cases 
in which the phases do not influence the final transmis- 
sion or reflection observed from a resonator, and therefore 
the coupling constants often are treated as real numbers, 
neglecting the phase information [3, [23| . 

Motivated by this, here, we study array of microres- 
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FIG. 1. (Color online) The considered setup of three coupled 
resonators in loop configuration. The resonators are probed 
by a fiber coupled to one of the resonators. Due to the ar- 
rangement of the resonators, light can evolve in a roundtrip 
process, e.g., from cavity 1 via cavities 2 and 3 back to cav- 
ity 1 without leaving the resonator array, which leads to rich 
interference effects. 
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onators in situations in which the phase of the coupling 
constants are crucial. In particular, we focus on situ- 
ations in which processes analogous to the closed-loop 
pathways in atoms occur. For this, we analyze an array 
of three coupled microresonators probed by a tapered 
fiber, see Fig. [TJ This setup is the simplest arrange- 
ment which allows for a closed-loop roundtrip pathway in 
analogy to the closed-loop atomic level systems. Due to 
the arrangement of the resonators, light can evolve in a 
roundtrip process, e.g., from cavity 1 via cavities 2 and 3 
back to cavity 1 without leaving the resonator array. We 
identify the evolution pathways for the photons entering 
the resonator array contributing to the transmission and 
reflection, and determine conditions for the dependence 
of these optical properties on the phases of the coupling 
and scattering constants. Based on these results, we an- 
alyze the phase dependence for several configurations in 
detail. Finally, we discuss possible applications for the 
phase-sensitivity. Our analytical results and interpreta- 
tions are based on quantum mechanical coupled mode 
theory. Additionally, we verify mechanism of the pro- 
posed applications using numerical finite-difference time- 
domain solutions of Maxwell's equations in two dimen- 
sion on a grid. 

This article is organized as follows. In Sec. [TT1 we de- 
scribe our model and the observables. In Sec. IIII1 we 
present our results. We start by analyzing the differ- 
ent pathways light can take through the resonator array, 
and how phase-sensitivity can arise from these pathways. 
We then move on to a discussion of transmission and re- 
flection spectra for phase-dependent systems, and of the 
most important roundtrip process enabled by the loop 
structure of the resonators. Finally, possible applica- 
tions are discussed in Sec. [IVJ Here, we show how to 
use our setup in order to measure the refractive index of 
a medium or to determine the position of a nanoparticle. 



II. THEORETICAL CONSIDERATIONS 



&i,im which is coupled into cavity 1. The scattering in- 
side cavity number n we describe by the parameter h n 
with n G {1,2,3}. The coupling between resonators n 
and m is described by the coupling constant £ mn with 
to, n G {1,2,3}. k denotes the coupling strength between 
the fiber and cavity 1. Using these definitions, in a suit- 
able interaction picture the Hamiltonian of our system 
reads (H 



H = H + Hl+ Hcs , 



with 



H = h^Ai(a]ai + b]bi) , 
i=i 

Hl = ihV2K[a ltln a\ - a* l in ai] , 
3 

HCS = ft J] + C m b \n a n) 



(1) 

(2a) 
(2b) 
(2c) 



Here for notational simplicity we defined h n — 
and used £ nm = £ mn . The detunings are defined as 
A; = uji — uji n , where w; is the resonance frequency 
of resonator I and u>i n is the frequency of the probing 
light. In the following, we assume equal resonance fre- 
quencies of the resonators Ai = A2 = A3 — A in our 
calculations. We write the complex couplings constants 



as 



I6» 



s l * m ™. For our calculations we assume 

for 



the critical coupling condition n = + (71 /2) 2 

the coupling between the fiber and cavity 1 to be ful- 
filled [40(. Since the modes {a m , b m } are assumed to have 
the same frequency it is reasonable to assume equal inter- 
nal loss rates 71 = 72 = 73 = 7 for all modes {a m , b m }. 
Then the total decay of mode a\ [&i] can be calculated 
according to j ai = 2k + 71 [7^ = 2k + 71], whereas 

7a 2 = lb 2 = 7a 3 = 7fc 3 = 7- 

From the Heisenberg equation, the time evolutions of 
the six mode operators a m and b m can be obtained as 



A. Description of the model system 

The system we consider consists of three equal nearby 
whispering gallery mode microresonators coupled to a ta- 
pered glass fiber, arranged as shown in Fig. [1] In such 
a toroidal microresonator a photon can be many times 
totally reflected at the edges of the cavity and thus move 
on a polygonal path very similar to a circle. The origi- 
nally undisturbed resonances of the cavity occur in pairs 
of clock- and anticlockwise propagating modes. Photons 
belonging to such a pair {a m ,b m } (m G {1,2,3}) differ 
only by their propagation direction but have the same fre- 
quency. Due to scattering processes, e.g. caused by ma- 
terial imperfections, the modes of such pairs can be scat- 
tered into each other. This scattering affects the eigen- 
modes of the system and thus also changes its eigenener- 
gies. 

As input field we consider a weak probe field of mode 



iA m + --f am ] a, t 



b m = 



H~ ^ml v^^^m^in 1 

1 



(3a) 



5A. 



;7fc„ 



bm-i^ Cm a n , (3b) 



where Sij is the Kronecker Delta function. Since only 
cavity one couples to the fiber, in the following we write 

O-l.in = O'ini a l,out = O'out an d fol.ojit = b ou f. 



B. Observables 

In our numerical calculations we investigate the trans- 
mission and reflection of light sent into the system via 
the coupled fiber. In particular, we are interested in the 
steady state mean values of the output mode operators. 
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We neglect fluctuations of the photon mode operators 
and calculate the steady state by setting a m = b m = 
in Eq. 03 For the calculation of the output operators, we 
use the input-output relation [4lj 



'2n(a. n 



\®out ) — Q>in 

(bout) =V2K(b m ) ■ 
Then the transmission and reflection become 
\(alut a out)\ 2 _ |(a «t)| 2 



T 

R = 



\Q"in | | 

\(bl ut bout)\ 2 ^ \(b out )\ 2 



(4a) 
(4b) 



(5a) 
(5b) 



III. RESULTS 
A. Light pathway analysis 

Light entering our three coupled resonator system can 
travel along a variety of different pathways before leav- 
ing cavity 1 as a transmitted or reflected photon. The 
different pathways have relative phases, since every time 
a photon is scattered inside one cavity or moves between 
two cavities the pathway amplitude is multiplied by the 
respective complex coupling constant or scattering pa- 
rameter with phase angle cj>ij. The total reflected and 
transmitted light then arises as the interference of all pos- 
sible pathways amplitudes. Therefore the phase angles 
4>ij can be expected to have a great impact on the result- 
ing output fluxes. In the following, we analyze conditions 
on the system parameters for obtaining phase-dependent 
or phase-independent reflections and transmissions. 

We start with a reduced system in which only phase in- 
dependent pathways are possible and afterwards explain 
how a phase dependence can arise. For this we set some 
of the couplings and scattering rates to zero. If we decou- 
ple cavities 2 and 3 from cavity 1 by setting £12 = £13 = 
as shown in Fig. UK A), no coupling process can take place 

but only scattering a\ — ^» b\ and vice versa is possible. 
Although the scattering parameter £n has an impact on 
the light's phase, the overall transmission or reflection 
are not phase dependent. The reason for this is that all 
pathways contributing to the transmission include no or 
an even number of scattering processes, e.g. from a\ to 
61 and back. Since the respective scattering parameters 
are complex conjugates, they have no net influence on 
the phase of the outcoming flux of each pathways. Simi- 
larly, the reflected light constitutes of contributions from 
pathways with an uneven number of scattering processes. 
This also means that all interfering contributions have 
the same phase when leaving the system and thus their 
interference is constructive independent of the phase of 
the scattering parameter. 

In the next step, we consider two coupled microcav- 
ities by setting £i 2 ^ whereas the scattering rate £22 





FIG. 2. (Color online) Example pathways of photons passing 
through the resonator array. In (A), the photon enters only 
the first cavity, scatters inside cavity 1, and then leaves the 
cavity back into the fiber. In (B), the photon traverses two 
cavities, but scatters only in cavity 1. (C) shows a pathway 
through two cavities with scattering in both the cavities, (D) 
a path through all cavities with scattering in cavity 1 only, 
(E) a path through all cavities without scattering, and (F) an 
evolution through three cavities with scattering in cavities 2 
and 3. 



in cavity 2 is assumed to be zero, see Fig. [U(B). Taking 
into account only small numbers of scattering or coupling 
processes, the following three different pathways are pos- 
sible: 

(6a) 
(6b) 
(6c) 

Note that we have omitted here pathways differing from 



ai 


€11 






ai 








a 1 


e 12) 


b 2 H 


ai — > bi 



4 



the ones in Eqs. © by a double scattering within cavity 
1 . In this reduced system, light evolving into cavity 2 al- 
ways returns to cavity 1 on the same way, and these pro- 
cesses are described by £12 and £* 2 , such that the phase 
of the coupling constant £12 does not affect the output in- 
tensities. Similar to the case in Fig. [21 A), also the phase 
of £n does not affect the transmission or reflection. 

Next, in addition we set £22 non-zero, see Fig. [^f C) . 
In this case, the phases cfiij do affect the transmission 
and reflection, since now in addition to the pathways in 
Eqs. ([5]), two additional leading order pathways are pos- 
sible: 
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Thus incoming light can propagate through cavity 1, en- 
ter cavity 2, scatter there, and then leave the system in 
reflection direction. Also, it can pass through cavities 1 
and 2 and be scattered in both cavities, leaving the sys- 
tem in transmission direction. These pathways depend 
on the phase of £22, £11 and £12, and thus the transmis- 
sion and reflection become dependent on these phases. 

Adding the third cavity, a qualitatively different evo- 
lution through the cavity system becomes possible. We 
first only chose the scattering rate £n as nonzero, see 
Fig. [2lD). Now incoming light can evolve on a roundtrip 
through all three cavities, be scattered in cavity 1 from 
mode b\ to a\ and afterwards leave the coupled resonator 
system in transmission direction. Processes like this in 
which light moves along pathways which describe a loop 
through all three cavities will play a special role in our 
further investigations, as they depend on the phases of 
all coupling constants between the coupled cavities. In 
the following we will refer to them as roundtrip processes. 
Beside these roundtrip pathways, also non-roundtrip pro- 
cesses are possible, and the transmission and reflection 
in leading order are superpositions of the pathways men- 
tioned in Eqs. © and the following pathways: 
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plus the pathway in which light enters cavity 1 and leaves 
the system without any further scattering. The resulting 
superposition of these amplitudes renders the transmis- 
sion and reflection dependent on both the phases of the 
couplings constants cfiij and the phases of the scattering 
rate (fin. 

It should be noted, however, that without scattering 
£ii = as in Fig. HfE) , also the dependence of the output 



intensities on the phases of the coupling constants cfiij 
disappears. All leading order pathways contributing to 
the reflection lead the light on a roundtrip pathway [see 
Eq. (|5a|) and (|8bl) ] such that they have the same final 
phase and thus interfere constructively for all choices of 
the phase angles cfiij . Similarly, the phase dependence in 
transmission direction vanishes. 

In the most general case with all coupling constants 
and scattering rates different from zero as indicated in 
Fig. |2JF), many interfering phase dependent pathways 
become possible, such that a dependence on all phases 
can be expected. In the following, we will analyze this 
phase dependence in detail. 



B. Transmission and reflection without roundtrip 
process 

We now turn to numerical results for the transmis- 
sion and reflection in the loop system. In Fig. [3] the 
transmission and reflection are shown for different values 
of the phase angles cfiij in dependence of the detuning 
A of the modes in cavity 1 to the incident light. We 
choose the scattering rates £n = 3O7, £22 = 20e i< ^ 22 7 
and £33 = 2O7. The couplings between the two cavities 
are £12 = 30e* 2 7, £13 = 30e* 3 7 and £23 = 0. Since 
£23 = 0, no roundtrip process is possible. For this choice 
of parameters, the pathways in Eq. ([B]) and ([7]) are possi- 
ble as shown in Fig.[HC), and additionally the analogous 
ones for cavity 3 instead of cavity 2. 

While studying the phase dependence, our variables 
are the angles </>i2, (fits and 022- In the solid lines in 
Fig. |3jA) all <fiij — 0, i.e. all coupling and scattering 
constants are taken as real numbers. In order to explain 
this result we consider the occupancy and the phase of 
each light mode inside cavity 1. Since coupling light out 
of the cavity into the glass fiber leads to a phase shift 
of 7T, the outcoupled light of mode a\ interferes construc- 
tively with aj n when its phase <fi ai = tt whereas we obtain 
destructive interference for cfi ai = 0, see also the input- 
output relations Eq. (|4]). As we choose £ R only 
the absolute value of <fi ai is of relevance for analyzing the 
interference. Therefore we define cfi a — \<fi ai \- Accord- 
ing to the input-output relations Eq. Q we can expect 
zero transmission T = if both the phase has a value 
leading to destructive interference in forward direction, 
and the amplitude satisfies 2k|oi| 2 = |eti„| 2 . If only one 
of the two conditions is fulfilled, only partial transmis- 
sion can be expected. In contrast, for mode b, there is 
no input field which can interfere with the field leaking 
out of cavity 1 into reflection direction. Therefore the 
reflected intensity is proportional to the intensity of b\ 
inside cavity 1. 

To verify this interpretation quantitatively, in the up- 
per two subfigures of Fig. [4] we show (A) the scaled pho- 
ton mode occupancies O ai and Ob 1 inside cavity 1 pro- 
portional to I ai 1 2 and |6i| 2 , respectively, and (B) the an- 
gle (fi a . Note that the amplitude condition for maximum 
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FIG. 3. (Color online) Transmission T and reflection R in the 
loop system. The parameters are £12 = £13 = 307 and £23 = 0, 
such that scattering between cavities 1 <-¥ 2 and 1 f> 3 are 
possible, but not between 2 and 3. Scattering occurs in all 
cavities with rates £11 = 3O7 and £22 = £33 = 2O7. In (A) 
all phases are chosen zero, <j>ij = 0. The dashed lines show 
corresponding results for transmission (T) and reflection (R) 
averaged over the phase angles $12 or 022- In (B), the phases 
are chosen as <^i2 = 0.2tt, and and all other 4>ij = 0. In (C), 



ii2 = — 0.6-7T and cf>i3 = QAn, and all other 



0. 



transmission corresponds to the horizontal line at scaled 
occupancies equal to one in Fig.[4jA). As one can see from 
FigBJA), the occupancies of both the modes a\ and 61 
have maxima around A = ±207. By contrast, the trans- 
mission in Fig. [3JA) shows minima for these detunings. 
The reason for this is that a (A = ±207) = 0. Thus, the 
field leaking from a\ into the fiber interferes destructively 
with the input field di n . Almost perfect suppression in 
forward direction T ss is achieved since for this detun- 
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FIG. 4. (Color online) Field configuration in resonator 1. 
(A) and (C): Scaled occupancies O ai = 2/-c|ai| 2 /|(2i n | 2 and 
Ob 1 = 2«;|6i| 2 /|ai n | 2 of the two counter-propagating modes 
ai and 61 inside cavity 1. (B) and (D): Phase cf> a = \4> ai \ of 
the counter-clockwise propagating mode inside cavity 1 which 
is coupled out into the fiber in transmission direction. In (A) 
and (B) , the parameters are chosen as in Fig. |3jA) , whereas 
in (C) and (D) they are as in Fig.E^B). 



ing 2«|ai| 



A similar interpretation holds for 



A = ±607. For resonant light, i.e. A = 0, a maximum in 
T can be observed in Fig. |3jA). This can be traced back 
to the low occupancy of mode ai in Fig 0|A) , such that 
the transmission mainly consists of a,i n . Note that for 
the parameters chosen in Fig. EJA), both the occupan- 
cies and the phases shown in Fig. are symmetric with 
respect to the detuning. From this it follows that also T 
and R are symmetric functions of A. 

In the next example in Fig.^B), we in contrast to (A) 
set </>i2 = 0.2-7T, but keep all other fyj = 0. Compared 
to the reflection and transmission in (A) , we can observe 
two additional peaks in the transmission and two addi- 
tional zeros in the reflection signal around A = ±207. In 
order to explain these results, we consider Fig. HJC) and 
(D). The minima in T at positions A = ±107, A = ±307 
and A = ±607 again arise from destructive interference 
with 4> a = and nearly fulfilled amplitude condition. At 
A = ±207 the modes of cavity 1 are nearly unpopulated. 
Thus R sa and the transmission is governed by the in- 
put flux cjj n . Interestingly, cavities 2 and 3 nevertheless 
contain much higher light intensities, which leads to de- 
coherence via 7. From T(A = ±20) = 0.82 ^ 1 we find 
that even though cavity 1 is almost empty, about 20% of 
the input light is coupled into the system. For larger de- 
tunings A all results in Fig. @|C) and (D) are similar to 
the graphs in (A) and (B). From this we conclude that in 
the offresonant case, as expected the phase angles <jf>y of 
the coupling and scattering parameters have only weak 
influence on the systems dynamics. Interestingly, the re- 
flection in Fig. E^B) is not a symmetric function of the 
detuning, while the transmission still is symmetric. This 
difference can arise since different pathways contribute 
to T and R. 

In the third example in Fig.^C) we chose ^12 = — 0.6-7T 
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FIG. 5. (Color online) Eigenvalue analysis to interpret the 
phase dependence. The left subfigures (A) and (C) show the 
imaginary part of the eigenvalues (£) of the matrix M de- 
scribing the system dynamics, which correspond to the en- 
ergy of the system's dressed states. They are plotted against 
the phase angle 4>i2- The right subfigures (B) and (D) show 
the power spectrum |£| 2 of one eigenvalue as example, reveal- 
ing the periodicity of the eigenenergies in 0i2. The parame- 
ters in (A) and (B) are fn = £13 = 3O7, £22 = £33 = 2O7, 
£12 = SOye 11 * 12 and £23 = 0, thus no roundtrip process is pos- 
sible. In (C) and (D) £23 = 157, thus the roundtrip process 
is possible. 

and 0i3 = 0.4-7T, but keep all other phases zero. Thus 
transitions between resonators lf>3 and 1 2 are pos- 
sible, but not between 2 and 3. It can be seen that this 
change in the phase of the coupling parameters leads to 
considerable modifications of the transmission and reflec- 
tion properties. Both transmission and reflection have a 
simple structure, but are not symmetric with respect to 
A. The interpretation of the peak structure is similar to 
the two previous cases. 

Analyzing the phase-dependence, we found that the 
transmission and reflection are 7r-periodic in <fii2 and </>i3 
but 2-7r-periodic in the phases of the scattering parame- 
ters 4>a. In the following section, we study and interpret 
this periodicity of our results in more detail using an 
eigenvalue analysis. 



complex eigenenergies of these dressed states. The real 
part of the eigenvalues of M can be interpreted as the 
decay rates whereas the imaginary parts correspond to 
the eigenenergies. Coupling constants as well as scatter- 
ing parameters included in M. shift the original eigenfrc- 
quencies of the resonators. Thus the imaginary parts £ 
of the eigenvalues of matrix M. are possible positions for 
peaks or dips in the reflection and transmission. These 
eigenvalues thus allow for a study of the dependence of 
the transmission and reflection spectra on the phase an- 
gles 4>ij. For example, to study the dependence on <p 12 , a 
Fourier transformation of one of the eigenvalues £ gives 



COM = 




The Fourier coefficients £(Z) then determine the peri- 
odicity of the eigenvalues in the phase 4>i2- First, we 
consider the case where no roundtrip process is possi- 
ble, i.e. £23 = 0. In Fig. [5] we show the six eigenvalues 
(A) and their power spectrum (B) for the parameters 
£11 = £13 = 307, £22 = £33 = 2O7, £12 = 3O7exp[z0i 2 ] 
an d £23 = 0. In this case, we observe that ((I) is different 
from zero only for even numbers of I. This means that the 
eigenvalues are 7r periodic in (f>i2 . If no roundtrip process 
can take place, light evolving from cavity 1 into cavity 2 
or 3 has to move the same way back in order to leave the 
system and to be detected as output light. Thus it inter- 
acts two times with the same coupling constant. This is 
the reason why the Fourier coefficient belonging to <pi2 
vanishes for uneven numbers of I. In Fig. [HC) and (D) 
we show the eigenvalues for £23 = 157- All further pa- 
rameters are the same as before. This setup corresponds 
to Fig. Op 1 ). In this roundtrip process is possible, 

and the Fourier coefficients £(7) are different from zero 
both for even or uneven numbers of /. The reason is that 
light can evolve, e.g., from cavity 1 via cavities 2 and 
3 back to 1 such that the phase (j>± 2 influences the path 
amplitude only once. Thus in this case, the eigenenergies 
are 27r periodic in (f>i2. 

Similar analysis allows to also reveal the phase depen- 
dence of the other coupling constants. 



C. Eigenvalue analysis 

In this section, we analyze the phase-dependence of 
the transmission and reflection spectra in a more general 
way. For this, we consider the eigenvalues of the matrix 
governing our system's dynamics. The equations of mo- 
tion for the six mode operators {aj, b{\ can be written as 

^-C = M ■ C , (9a) 
C = (01, h, a 2 ,b 2 ,a 3 , b 3 ) T . (9b) 

By diagonalizing M. , the dressed states of the system can 
be evaluated. The complex eigenvalues correspond to the 



D. The roundtrip process 

In this Section we study the roundtrip process and the 
resulting effects in detail. For this, we make use of the 
fact that by turning the coupling £23 on and off it can 
be controlled whether or not a roundtrip pathway can 
be taken by entering light. Therefore, in the first step 
we expand the transmission T in this coupling constant 
around £23 = 0. This expansion will reveal the leading 
order effects of the roundtrip process for small £23- Af- 
terwards we present numerical results for T and R for 
more general parameter sets which allow the roundtrip 
processes to take place. 
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FIG. 6. (Color online) Expansion of the transmission in or- 
ders of the roundtrip process. For this, the transmission is 
expanded around in |^23 1 /t- The parameters are £11 = 5O7, 
£22 = 20e lte 7 , ^3 = 20 7 , £23 = IO7, fis = 30e^ 13 7 and 
all <f>ij = 0. (A) zero order contribution (£23 = 0), (B) first 
order contribution ci, (C) second order contribution C2, (D) 
complete transmission for £23 = 37. 



1. Expansion in orders of the roundtrip process 

To elucidate the impact of the roundtrip process on 
the transmission spectrum, we performed a Taylor expan- 
sion of the transmission amplitude to the second order in 
I £23 1/7 around £23 = 0. This refers to the case where the 
coupling between cavity 2 and 3 is much weaker than the 
other couplings The expanded transmission reads: 



|(q nt)| s 

\(a in )\ 2 



CO + Cl 



16 



7 



c 2 - 



I63I 2 



r 



(11) 



Here, c„ are the Taylor expansion coefficients. The 0th 
order corresponds to the case where no roundtrip process 
is possible, i.e. £23 = 0. The respective result is shown 
in Fig. [5] (A). The parameters are chosen as £n = 5O7, 
62 - 20e* 2 7, 63 = 20 7 , £12 = 10 7 , = 30e**"7 
and all <f>ij = 0. In Fig. [6jB) and (C) we show the first 
and second order Taylor coefficients c\ and C2 which pro- 
vide the respective correction terms for the transmission. 
The first order corresponds to a pathway in which light 
interacts once with the coupling constant £23- The most 
probable process of this kind that ends in a\ is a single 
roundtrip through all three cavities, see e.g. Eq. (|5c| 
and (18d[) . The second order contributions correspond to 
pathways where light moves twice between cavity 2 and 
3, including a double roundtrip process. Examples for 
such pathways are given in Eq. (|5e|) and (|8f|) . 

We can see from these figures that the first order cor- 
rection leads to a small dip that overlaps with the broad 
resonance round A = 0. By contrast, the second or- 
der correction consists of two sharp resonances round 
A = 2O7 = £22 = £33- Comparing the sum of the curves 
of Fig. [11(B) and (C) to the full output for $23 = 37 (D), 




FIG. 7. (Color online) Transmission (blue) and reflection 
(red) spectra for parameters which allow for the roundtrip 
process. The parameters are the same as in Fig. [6] except for 
£23 = 57. In (A) and (B), all phases are zero (cj>ij — 0). In 
(C) and (D), 022 = 1.6-7T and ^13 = OAn and all other phases 
are chosen zero. In (E) and (F), ^22 = 1-67T and <j>i3 — — OAn 
and all other <j>a =0. 



we find that for the small value of £23 taken in our ex- 
ample the first two orders in the expansion are sufficient 
to almost perfectly approximate the full result in (D). 

The structure of T can be explained using the same 
interpretation techniques as applied for Fig. [3] based on 
the occupancy and phase of the fields inside cavity 1. 



2. Transmission and reflection with roundtrip process 

Next we study the phase dependence of the trans- 
mission T and reflection R for a choice of parameters 
for which a roundtrip process is possible. For this, we 
choose all coupling constants and all scattering parame- 
ters non-zero, as shown in Fig. [2JF). All parameters are 
as in Fig. [5] except for £23 = 57. We start by consid- 
ering the transmission and reflection if all phase angles 
4>ij — 0. The respective results are shown in Fig. EJA) 
and (B). We can see clearly the two side band dips around 
|A| = 2O7 = |£22 1 = I £33 1 arising from the second order 
Taylor correction of Eq. (fTTj) . 

Fig. E|JC) and (D) show corresponding results with 
phase angles changed to ^22 = l-fwr and ^13 = OAn. We 
observe that the reflection becomes asymmetric and only 
one sharp dip around A = 2O7 remains. The transmis- 
sion also changes, but remains symmetric. Upon chang- 




FIG. 8. (Color online) Transmission (blue) and reflection 
(red) spectra for parameters which allow for the roundtrip 
process. The spectra are plotted against the phase (j>i3- The 
other parameters are chosen as in Fig. [7] with A = —19. 57. 



ing the sign of phase 0i3, the sharp dip moves to the 
opposite side of the spectrum, i.e., to A w — 2O7, see 
subfigures (E) and (F). The transmission still remains 
symmetric, but is also affected by the phase change of 
013. In Fig. [8] we show T and R in dependence on the 
phase angle 4>\3 for A = —19.57 which is the position of 
one of the narrow structures. We find that even for small 
changes of the angle (^13, the transmission and reflection 
can changes considerably. 

We thus conclude that the possibility of taking a 
roundtrip pathway in the loop system is the origin of 
narrow structures in both the transmission and the re- 
flection, and these narrow structures are sensitive to the 
phases of the coupling constants. This invites applica- 
tions based on the dependence of these coupling con- 
stants on an observable. In the following Section we de- 
scribe two possible applications based on the sensitivity 
of the transmission spectrum on the phase angles ^ij ■ 



IV. APPLICATIONS 

In this Section we discuss two possible applications of 
our setup. They rely on the dependency of the trans- 
mission and reflection on the scattering and coupling pa- 
rameters. First, we analyze the possibility to detect the 
position of a small particle such as an atom or a nano ob- 
ject close to one of the resonators as shown in Fig. [^A). 
Second, we aim at measuring small changes in the refrac- 
tive index between two cavities. These changes could be 
induced by an object placed in the free space, or by em- 
bedding the interface area in a liquid. The corresponding 
setup is sketched in Fig. 

In the following, we first provide a theoretical back- 
ground to our calculations, and then discuss numerical 
results both based on the coupled mode theory and on a 
numerical integration of Maxwell's equations on a grid. 
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FIG. 9. (Color online) Setups for the two example applica- 
tions. (A) shows the measurement of the position of a particle 
at position R p close to resonator 2. (B) shows the measure- 
ment of changes in the index of refraction between resonators 
2 and 3 induced, e.g., by slabs of different refractive indices. 

A. Monitoring the position of a nano particle 

A subwavelength refractive object located very close 
to one of our cavities as shown in Fig. A) gives rise to 
a scattering of the fields propagating inside the cavities. 
This results in damping and in a c oupling between conter- 
propagating WGM pairs {a n , b n } [42|, |43|. The scattering 
resulting from a nano particle has been studied in [42l . l44l |. 
where it was demonstrated that the size of a nano par- 
ticle can be determined using high-Q WGM. This raises 
the question, whether also the particle position can be 
determined. To address this question, we note that in 
the coupled mode theory, the complex electric fields of 
the modes a n and b n in the n-th cavity are related and 
can be written in cylindrical coordinates R = (p, 9, z) T 
as [3 

Eg(R) = (£>, z),iE$(fi, z),E° z (p, z)f e™ e (12a) 

E°(R) = (£>, 2), -iE° e (p, z),E° z (p, z)f e~ m9 . 

(12b) 

The origin of the used coordinate system is at the cen- 
ter of the WGM cavity. For high-Q WGMs with small 
loss rate, a n and b n are to a good approximation com- 
plex conjugates of each other, such that the three com- 
ponents Ep,Eg,E z are real functions (44J. We assume 
that a scattering particle at position R p giving rise to a 
point-like dielectric fluctuation with size much smaller 
than the wavelength. We denote the dielectric con- 
stant of the resonators as e c , that of the homogeneous 
medium surrounding the resonators as s s , and that of 
the medium with added particle as e p (R). The scatter- 
ing parameter £ nn is the proportional to the difference 
fe(R) = £ p (R) — e s as well as to the intensity of the 
electric field at the position of the particle, and can be 
written as @, Hj| 

_ u m Jy p MR) - £ s )E°*(R)Eg(R)rfR 
Un ~ 2 / £s |E£(R)PdR (13a) 

oc fc(Rp)e 2lm9 " (R p ) |E°(R p )| 2 . (13b) 

Here, m is the azimuthal mode number, and in Eq. (|13bl) . 
for simplicity, we set E°(R) = Eg(R) = E°*(R). It can 
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be seen that the scattering resulting from the particle 
is a complex number. Its phase depends on the posi- 
tion n (R p ) of the particle, see Fig. EIA). According to 
Eq. (|13b[) . this angle enters the phase of the coupling 
constants via 2m(9 n (R p ), such that an increase of the az- 
imuthal mode number leads to higher position sensitivity, 
but at the cost of a smaller range of uniquely determined 
positions, since the phase is only determined modulo 2tt. 
Since we found in the previous sections that the trans- 
mission and the reflection in our loop setup is sensitive 
to the phase of the coupling constants, in principle, a 
position determination becomes possible. 



B. Monitoring the dielectric constant of a nano 
slab 



Next, we turn to the measurement of the dielectric 
constant of a thin object located in the space between 
two of the cavities. Alternatively, this slight change could 
also be induced by the concentration of a fluid between 
the two resonators, or by varying the temperature [46| . 
In general, modifying the dielectric constant in between 
the two cavities gives rise to a change in the coupling 
of two cavities and to scattering. However, for a larger 
sample exceeding the wavelength scale, the scattering can 
be small such that the change in the coupling constant is 
dominant. In the following, we assume this condition to 
be fulfilled and neglect the scattering induced by the slab, 
and consider a slab in the region V s i a b- The total coupling 
can then be separated into three parts: The coupling 
without slab £im, the contribution £nm^ from the slab 
with a reference constant e° slab inserted in the gap of two 

cavities, and a change S^nm^ as the dielectric constant 
of the slab varies according to 6e s i a b = £siab~ e 1i a b- Thus 
the total coupling is given by 

$l=M [ (e c - £s )E°*(R)E^(R)dR (14a) 

JV cavU y 

e n t ab) =M [ (e° slab - £s )E°*(R)E^(R) dR (14b) 

JVslat 

= W (esiab - £ L b )E°*(R)E^(R)dR 

Jv„„ h 



t — <K°) I Aslab) , Si Aslab) 
snm Swn ~ Snm ~ "Snm > 

where 



N= ^T\j £ *l E °( R )l 2dR / £,|E^(R)| 2 rfR 



(14c) 
(14d) 

-1/2 

(15) 



Using a similar approximation as in case of the nano par- 
ticle, the change djnm^ induced is proportional to Se s i a b 
and given by [4J, |45j 



(slab) 
nm 



£ 



slab t(slab) 

£— ?nm 

lab fc s 



(16) 



Thus a large static coupling £nm is favorable. Again, 
since the transmission and reflection in our setup depend 
on the coupling constants, such a variation of the cou- 
pling constants can be detected. 



C. Reflection and transmission averaged over 
coupling and scattering phases 

In this Section we study the impact of an averaging 
over certain phase angles 0jj on the transmission and re- 
flection spectra. In all cases, we average over the full 
range of 2tt. The averaging over the coupling constant 
phase 0i2 could be visualized as an experimental setting 
in which the refractive index of the medium between cav- 
ity 1 and 2 changes between several measurements, e.g., 
due to changes in the concentration of a fluid filling this 
region. The averaging over a scattering constant 022 can 
be visualized as arising from different particle positions 
throughout the measurements [47] . 

In general, we find that the averaged curves differ con- 
siderably from the curves obtained for fixed phases such 
as (f>ij = 0. Fig. EIA) shows that averaging over i2 or 
4>22 lead to the same results, as long as resonators 2 and 
3 are uncoupled, i.e., £23 = 0. This can be understood 
by noting that in case of £23 = 0, all pathways in which 
012 or 2 2 lead to a final phase shift of the corresponding 
amplitude include both an interaction with the coupling 
£12 and an interaction with the scattering rate £22, see 
Eq. ([7a|). Therefore either averaging has the same effect. 

In contrast, if £23 ^ 0, the spectra averaged over the 
scattering phase 022 or over the coupling phase 0i2 dif- 
fer, see Fig. [TUJ The reason is that now 0i2 and 022 
affect different pathways. For example, the pathways in 
Eqs. (f8aj) -(|8d j) include a phase contribution of 0i2, but 
not of 022- Thus in contrast to the case without roundtrip 
process, the averaging over the two phases leads to dif- 
ferent results. 



D. Numerical simulations with FDTD 

In this section we numerically verify the suggested ap- 
plications using a finite-difference time-domain (FDTD) 
solution of Maxwell's equations on a grid [48]. The sim- 
ulations consider two dimensions for TM modes propa- 
gating in x direction (field components E x , E y , H z ) on a 
Yee-grid with grid size of 30 nm. The grid boundaries are 
modeled using Berenger-type perfectly absorbing bound- 
ary conditions, and we use a point-like sinusoidally os- 
cillating source in H z . We consider a background with 
permittivity e s = 1 , and resonators and a waveguide with 
permittivity e c — 4. The waveguide resonators have 
outer radius 3500 nm and inner radius 3350 nm. The 
distance between waveguide and resonator is 120 nm, 
the distances between the resonators are chosen equal 
as 200 nm. The waveguide has a width of 150 nm. 
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FIG. 10. (Color online) Transmission T, T and reflection _R, R 
for parameters as in Fig. [5] In (A), the solid lines T, R show 
the case with phases <j>ij = 0, the dashed lines T, R show 
results averaged over 022 ■ In (B), the solid lines T, R show 
the case with phases (j>ij = 0, the dashed lines T, R show 
results averaged over 0i2. 



A typical example for three resonators in loop config- 
uration is shown in Fig 111! The figure shows the field 
component H z after the time evolution has reached a 
stationary state. The pointlike source is in the lower left 
corner, and is placed in the center of the waveguide which 
runs along the lower edge of the figure. Since the point- 
like source does not exclusively excite waveguide modes, 
circularly spreading background waves originating from 
the source can be seen as well. The excitation then runs 
along the waveguide to the three resonators seen as the 
circular field arrangements. In this particular example, 
the two upper cavities contain standing wave excitations, 
which manifest themselves a modulated total intensity 
( "blinking" ) in the time-dependent dynamics of the field 
configuration. The right half of the lower cavity exhibits 
a less pronounced standing wave, whereas the left half 
is mostly filled by a running wave in clockwise direction. 
This "blinking" can be seen by comparing the figure to 
a corresponding snapshot slightly later in time. If the 
time is chosen appropriately, the time evolution of the 
standing waves is close to a minimum, such that the 
bright field regions in Fig Qj] are almost invisible. In con- 



FIG. 11. (Color online) Configuration of the field component 
H z obtained from the FDTD-simulation for three resonators 
in loop configuration without particle or slab. The image 
shows one snapshot in time. The excitation wavelength is 
A = 571.8 nm. See text for further explanation. 
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FIG. 12. (Color online) Transmission spectrum obtained from 
FDTD simulations for three resonators in loop configuration 
without particle or slab. 

trast, the running wave parts remain similar. In total, 
the setup in Fig QT] leads to a weak forward transmission 
(T 1), which can be seen from the low field excita- 
tion downstream of the resonators. The energy is instead 
mostly reflected, which again is evidenced by standing 
wave field components in the waveguide between source 
and resonators. In contrast, the other waveguide parts 
only carry running wave excitations, as expected. 

To evaluate the transmission, we sum the time- 
averaged Poynting vector contributions in a plane 
transversal to the waveguide to the right of the res- 
onators. The obtained energy flux is normalized to the 
case without resonators. By varying the frequency of the 
excitation, a transmission spectrum can be determined. 
A typical example is shown in Fig. 1121 

For the applications, we add a particle or a slab to 
the cavities as indicated in Fig. [S] The particle has a 
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FIG. 13. (Color online) Transmission spectrum obtained from 
FDTD simulations for three resonators in loop configuration 
with a particle placed close to the top left resonator. The 
different curves show particle positions (i) 6 — 90°, (ii) 6 = 
95°, and (hi) 6 = 180°. 
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FIG. 14. (Color online) Transmission spectrum obtained from 
FDTD simulations for three resonators in loop configuration 
with a slab placed between the two rightmost cavities, see 
Fig. HfB). (A) and (B) show two spectral lines at different 
wavelengths, (i) corresponds to a slab with permittivity e = 
4.0, (ii) to a slab with e = 4.1. 



radius of 90 nm, and it is placed at different azimuthal 
positions with a radial distance of 90 nm to the top left 
cavity. The transmission spectrum around one spectral 
peak is shown for different particle positions in Fig 1131 
Other than the azimuthal position, no parameters are 
changed between the three curves in this figure. It can 
be seen that changing the position leads to a shift of the 
resonance line in the transmission spectrum. This shift 
can be observed over a large range of position values, as 
indicated by the lines with 90 degrees and 180 degrees. 
On the other hand, a shift can already be observed for 
a change of position by 5 degrees. This has to be com- 



pared to the angular range of about 7 degrees for one 
wavelength in the azimuthal mode number 52 realized in 
this numerical example. 

Fig. [TH shows corresponding results for a slab. The 
slab has width 60 nm and length 1500 nm and is placed 
symmetrically between the two rightmost resonators. We 
consider slabs with permittivity e = 4.0 and e = 4.1. 
The two subfigures in Fig. [TJ] show two resonances in the 
transmission spectrum. It can be seen that the change 
in the permittivity from 4.0 to 4.1 shifts the position of 
a resonance line. It is interesting to note that not the 
whole spectrum is shifted, but only part of the spectral 
lines, as can be seen by the parts in Fig. [14] which are 
not affected by the change of the slab. 

We thus conclude that the FDTD simulations of the 
considered loop systems with additional particle or slab 
serve as a proof or principle for the applications discussed 
in Sec. IIVI While a direct connection to the coupled- 
mode calculations is not possible since the individual cou- 
pling constants and resonator properties realized in the 
numerical simulations are unknown, both the variation 
of the particle position and the index of refraction of the 
slab led to a shift of resonance lines in the transmission 
spectrum, and therefore should be detectable. 

V. SUMMARY 

We analyzed interference effects in a system consisting 
of three coupled microcavities arranged in such a way 
that light can evolve through the resonators in a non- 
trivial loop roundtrip. The system is probed by a fiber 
coupled to one of the resonators. The interplay of the 
different pathways light can take while passing through 
the resonator array leads to rich structures in the trans- 
mission and reflection spectra of the system. In partic- 
ular, we have focused on a sensitivity of the spectra on 
the phases of the different scattering and coupling con- 
stants. We found that the roundtrip process in which 
light moves on a circle through all three cavities, enabled 
by the special arrangement of the resonators, leads to 
additional pathways which increase the sensitivity of the 
spectra to the phases. Finally, we discussed two appli- 
cations for the found phase-sensitivity. First, we studied 
the determination of the position of a particle placed in 
the evanescent field of one of the resonators. Second, 
we analyzed the measurement of the index of refraction 
of a slab placed between two cavities. Our results are 
on the one hand based on quantum mechanical coupled 
mode theory, which allows to interpret the spectra in de- 
tail based on the underlying physical mechanisms. On 
the other hand, we verified the sensitivity of the spectra 
in the two discussed simulations using numerical finite- 
difference time-domain simulations of Maxwell's equa- 
tions on a grid. 
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